Dynamics of Kinks: Nucleation, Diffusion and Annihilation 
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We investigate the nucleation, annihilation, and dynamics of kinks in a classical (l + l)-dimensional 
4> 4 field theory at finite temperature. From large scale Langevin simulations, we establish that the 
nucleation rate is proportional to the square of the equilibrium density of kinks. We identify two 
annihilation time scales: one due to kink-antikink pair recombination after nucleation, the other 
from non-recombinant annihilation. We introduce a mesoscopic model of diffusing kinks based on 
"paired" and "survivor" kinks/antikinks. Analytical predictions for the dynamical time scales, as 
well as the corresponding length scales, are in good agreement with the simulations. 
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PACS numbers: 05.20.-y, ll.10.-z, 63.75,+z, 64.60.Cn 



Many extended systems have localized coherent struc- 
tures that maintain their identity as they move, interact 
and are buffeted by local fluctuations. The statistical 
mechanics of these objects has diverse applications, e.g., 
in condensed matter physics [1], biology [2], and parti- 
cle physics [3]. The model to be studied here is a kink- 
bearing (f) 4 field theory in (1 + 1) dimensions, popular 
because its properties are representative of those found 
in many applications. Static equilibrium quantities of 
this theory, such as the kink density and spatial correla- 
tion functions, are now well understood and recent work 
has shown that theory and simulations are in good agree- 
ment [4-6]. However, dynamical processes, both close to 
and far out of equilibrium, are much less well understood. 
Questions include: What is the nucleation rate of kink- 
antikink pairs? How is an equilibrium population main- 
tained? How do these dynamical processes depend on 
the temperature and damping? These questions, among 
others, are the subject of this Letter. 

We introduce and analyze below a simple model of 
kink diffusion and annihilation that predicts the nucle- 
ation rate and provides a picture of the physical situa- 
tion, including the existence of multiple time and length 
scales. We also carry out high resolution numerical sim- 
ulations. As one consequence of our work, we are able to 
settle a recent controversy as to whether the nucleation 
rate of kinks in an overdamped system is proportional to 
exp(-2E k /3) [7] or exp(-3£ fc /3) [8] in favor of the first 
result (E k is the kink energy and (3 = l/k B T). 

We consider the dynamics of the 4 field obeying the 
following dimcnsionless Langevin equation [4] : 



d? t( f> = di x <t> + 0(1 - 2 ) - ^ + e(x,t). 



(i) 



with the fluctuation-dissipation relation enforced by 
(£(x,t)£(x',t')) = 2r)0- 1 6(x - x')6(t - t'). We perform 
simulations on lattices typically of 10 6 sites, using a finite 
difference algorithm that has second-order convergence to 
the continuum [6] . Typical values of the grid spacing and 
time step are Ax = 0.4 and At = 0.01. 



At zero temperature, the static kink solution centered 
at x = xo is <pk{x) = k(x — xo) where k{x) = tanh(x/v / 2); 
the corresponding antikink solution is 4> a {x) = —k(x — 
xo). Because there are only two potential minima, kinks 
alternate with antikinks on the spatial lattice. Imposing 
periodic boundary conditions constrains the number of 
kinks and antikinks to be equal. During the time evo- 
lution, we identify kinks and antikinks individually and 
follow the "lifeline" of each kink or antikink (Fig. 1). 



5.0 - 



4.5 - 




1 .0x1 



2x1 



FIG. 1. Timelines of kinks and antikinks: A small space- 
time portion of one numerical solution is shown. Many re- 
combinant nucleation-annihilation events are barely visible, 
forming small closed loops. /3 = 7, r] = 1. 

Equilibrium properties of one-dimensional systems, 
such as the free energy density and the correlation func- 
tion (4>(0)(f)(x)), can be calculated using the transfer in- 
tegral method [9]. The calculation is exact, although 
one typically must evaluate eigenvalues of the resulting 
Schrodinger equation numerically. When the on-site po- 
tential has the double-well form, as is the case here, one 
part of the free energy density at low temperature can 
be interpreted as due to kinks, forming a dilute gas with 
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density [9]: pu oc ^/E k {3exp(— E k (3). This WKB approx- 
imation is consistent with recent simulations at j3 > 6, 
where unambiguous identification of kinks is possible [4] . 
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FIG. 2. Nucleation rate, measured from numerical solution 
of (1) with tracking of timelines. The rate, F, of production 
per unit length per unit time of kink-antikink pairs is shown 
versus (3 for three values of r\: rj — 0.2 (triangles), r\ = 1 (filled 
circles) and rj = 5 (diamonds). The solid line is F = p\ and 
the dashed line is the best fit to F = ap\ for 77 = 1. 

An equilibrium density of kinks is maintained by a dy- 
namical balance of nucleation and annihilation of kink- 
antikink pairs (Fig. 1). The dependence of the nucle- 
ation rate Y on temperature and damping, however, is 
not directly calculable from the transfer integral; nor are 
unambiguous results for symmetric potentials available 
from saddle-point calculations [10]. While analogy with 
the Kramers' problem suggests Y oc exp(—2(3E k ) [11], the 
relationship T oc exp(— 3f3E k ) has also been suggested [8]. 
Our direct counting of nucleation events establishes that 
their rate is proportional to the square of the equilibrium 
density, that is T oc exp(— 2(3E k ) (Fig. 2). Below we show 
how this relation can be understood from a mesoscopic 
model of diffusing kinks with paired nucleation. 

At equilibrium, the nucleation rate is related to the 
mean kink lifetime r by p k — FY. Previous attempts to 
evaluate Y numerically [12,13] have proceeded by count- 
ing the number of kinks n k (t) and assuming exponential 
decay of (n k (t + T)n k (t)}. Unfortunately this approach 
provides no information on the underlying processes, and 
yields incorrect results if kinks are not properly identified 
on the lattice. In particular, results that appeared to sup- 
port r oc exp(— 2>(3E k ) were performed at temperatures 
too high for accurate computation of (n k (t+T)n k {t)) [12]. 

Because we identify individual nucleation events and 
follow individual kink lifelines, we can distinguish 
"paired" kinks (whose partner antikink is still alive) 
from "survivor" kinks (whose partner has been killed). 
We also distinguish and measure the contributions to 
the annihilation rate from the recombinant and various 
non-recombinant mechanisms (Fig. 3). The most fre- 



quent annihilation event is recombination of a recently- 
nucleated pair (designated I in Fig. 3) [7]. However, the 
"survivor" kinks that remain after a non-recombinant an- 
nihilation event (II or III) have a longer mean lifetime. 
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FIG. 3. The four kink-antikink annihilation processes. I: 
Recombination of a paired kink and antikink; II: Annihilation 
of a survivor with a paired kink/antikink; III Annihilation of 
a kink and an antikink from two neighboring nucleation sites 
(producing a survivor kink and antikink); IV Annihilation of 
two survivors. 

At finite temperature, the mean-squared displacement 
of an isolated kink is given by (Xj) = 2Dt. The diffu- 
sivity D can be estimated by using the zero-temperature 
kink as an ansatz in the equation of motion (1), yielding 
D ~ (EkPrf)- 1 [14], where E k = y^8/9 for a static kink. 
Corrections to D, arising because of fluctuations in the 
kink shape, are proportional to /3~ 2 and subdominant in 
the temperature range considered here. 

Our numerical observations, in particular that kink- 
antikink collisions at moderate to large damping always 
result in annihilation, motivate us to introduce the fol- 
lowing mesoscopic model of kink dynamics: (i) kink- 
antikink pairs are nucleated at random times and po- 
sitions with initial separation b <C pZ \ (ii) once born, 
kinks and antikinks diffuse independently with diffusiv- 
ity D; (iii) kinks and antikinks annihilate on collision. 
The separation between a kink and its partner performs 
Brownian motion with diffusivity 2D. Thus, if only re- 
combinate annihilation (I in Fig. 3) were allowed, the 
time to between nucleation and annihilation would have 
the density f t V[t <t]= bt~i (8irD)-i exp(-g^) [15]. 

To analyse our model, we use the following approxima- 
tion for non-recombinant annihilation: as long as both 
members of a pair are alive, there is a constant probabil- 
ity \x per unit time of a member being struck and "killed" 
by an outsider, i.e. of an event II or III. Thus, to each 
pair we assign a killing time , distributed according to 
"P[t M > t] = exp(-fit). Non-recombinant annihilation 
happens with probability Pft^ < t ] = 1 — exp(— bv), 
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where v 2 = [16,17]. The killing rate p depends on 
the density of kinks; we estimate it as follows. A new- 
born pair finds itself in a domain between an existing kink 
and antikink of typical length l/(2pk). The mean time 
for a diffusing particle to exit the region is proportional 
to (2D/?2)-i. Therefore, let 

f i = 2Da 2 pl (2) 

The value of the dimensionless factor was obtained from 
numerical measurements of length and timescales (see 
Figures 4 and 6 below): we estimate a ~ 8. 



Let R(t) = 4rV[t < t\t < tj. Then 



R(t) 



*(&)exp(-— )r 



exp(-/zi), 



(3) 



where N(b) = b(8irD)~2 exp(ub). In Fig. 4 we plot (3) 
and a histogram of values of t obtained from a large 
numerical solution of (1). The behavior R(t) oc t~5 is 
characteristic of Brownian excursions [15]. 
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FIG. 4. Recombination time density. Only recombinant 
histories (I in Fig. 3) are counted. The solid line is (3). The 
histogram is from a simulation of (1). /? = 6, 7? = 1. 



Our mesoscopic model has two timescales [17]: 
/ \ ^ 

TO = (t |t < t M ) = 
7> = (V|t M < t ) = 




(4) 
(5) 



The mean recombination time (4) depends on b; in con- 
trast r M has a non-zero limit for vb — » 0: t m — > l/(2/x). 
With the approximation that a "survivor" kink/antikink 
has the same probability per unit time, /x, of collision and 
death, the mean lifetime of a kink or antikink is given by 

1 



T e 



-vb 



-vb 



) + ^(l-c-). (6) 



As pb -» 0, r -» (3/2)6(2/iD)- 1 /2 = (3/4)6(aDp)- 1 . In 
the same limit, combining (2) and (6) yields the predic- 
tion that the nucleation rate is proportional to the square 
of the equilibrium density: 



Pk/r = ^ Da Pk, 



(7) 



The relation r oc p\ is also found in the discrete-space 
Ising model [18]. In contrast, the nucleation rate is pro- 
portional to p\ in systems where nucleation does not oc- 
cur in pairs [18,19]. The latter scaling was incorrectly 
predicted for the 4> 4 system [8], from an estimate of the 
annihilation rate that does not take into account paired 
nucleation. (In the </> 4 system one does, however, find 
that the rate of survivor-survivor annihilation events - 
IV in Fig. 3 - is proportional to p\.) In the </> 4 SPDE, 
the parameters D, T and b have in general a weak (non- 
exponential) temperature dependence. The lcngthscale b 
is of the same order as the width of an isolated kink. 
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FIG. 5. Nucleation rate versus damping for f3 = 6. The 
dashed line has slope I/77. 

In Fig. 5 we plot the measured nucleation rate versus 
the damping coefficient at fixed temperature. The nucle- 
ation rate is proportional to rj^ 1 for n 3> 1 [in agreement 
with (7)] and appears to plateau for 77 — > 0. At low damp- 
ing, however, direct measurement of the nucleation rate 
is problematic because kink-antikink collisions may re- 
sult in single or multiple bounces rather than immediate 
annihilation [20]. 

We now turn to the lengthscales in the system. A his- 
togram of distances between neighboring kinks and an- 
tikinks is well-approximated by an exponential with char- 
acteristic length (2pk)^ 1 . This simple form results from 
the cancellation of the tendency of paired kinks to be 
closer together than (2pk)^ 1 with the opposite tendency 
of survivor kinks. In Fig. 6 we plot f(x) =(number of 
occurrences of separation £ (x, x + dx))/(L dx). We also 
construct the histogram for the separations of only paired 
kinks and antikinks. The dashed curve is the probability 
of being at x, averaged over the lifetime, for a Brownian 
motion killed at x — and at rate p [16]: 



l(x) = Aexp(— vx) = N cxp(-apkx). 



(8) 
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FIG. 6. Kink-antikink spacing. In the upper part of the 
plot, a histogram of all next-neighbor spacings is plotted with 
(2pk) 2 exp(— 2pkx) (solid line). The lower part of the plot is 
a histogram of kink-antikink spacings only from paired kinks; 
the dashed line is (2p k ) 2 exp(—8pkx). (3 = 6, r\ = 1. 

The classification of kinks into paired kinks and sur- 
vivors, with the approximation that kinks have a con- 
stant probability (2) per unit time of non-recombinant 
annihilation, allows us to construct a macroscopic rate 
theory for the two densities n p (t) (paired kinks) and n s (t) 
(survivor kinks). The equation for n p (t) has a positive 
term due to nucleation and a negative term inversely pro- 
portional to the lifetime of pairs, (4). The terms in the 
equation for n s (t) correspond to processes III and IV in 
Fig. 3. Note that process II does not change the number 
of survivor kinks. We obtain 

hp = r — 2b~ 1 aD(n s + n p )n p 

h s = Da 2 (n s + n p ) 2 (n p - 2n s ). (9) 

The steady state solution of (9) gives the relationship (7) 
between Y and the equilibrium kink density. Nonequilib- 
rium dynamics are also correctly described: if Y = 0, the 
paired density quickly decays and, for late times, h s oc 
n^, in agreement with an exact result for the survival 
probability in the diffusion- limited reaction A+A^O [21]. 
While not exact, (9) illustrate that at least two coupled 
equations are necessary to capture the two timescales in 
the dynamics: no single rate equation can suffice. 

We have benefited from discussions with Kalvis Jan- 
sons, Eli Ben-Nairn and Vincent Hakim. Computations 
were performed at the National Energy Research Sci- 
entific Computing Center (NERSC), Lawrence Berkeley 
National Laboratory. 
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